3.31.63 \(\int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx\) [3063]

3.31.63.1 Optimal result
3.31.63.2 Mathematica [A] (verified)
3.31.63.3 Rubi [A] (warning: unable to verify)
3.31.63.4 Maple [A] (verified)
3.31.63.5 Fricas [F(-1)]
3.31.63.6 Sympy [F]
3.31.63.7 Maxima [F]
3.31.63.8 Giac [A] (verification not implemented)
3.31.63.9 Mupad [F(-1)]

3.31.63.1 Optimal result

Integrand size = 22, antiderivative size = 135 \[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=-\frac {3 b d \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}}{2 a^2 \sqrt {\frac {d}{x}}}+\frac {\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}} x}{a}-\frac {\left (4 a c-3 b^2 d\right ) \text {arctanh}\left (\frac {2 a+b \sqrt {\frac {d}{x}}}{2 \sqrt {a} \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}}\right )}{4 a^{5/2}} \]

output
-1/4*(-3*b^2*d+4*a*c)*arctanh(1/2*(2*a+b*(d/x)^(1/2))/a^(1/2)/(a+c/x+b*(d/ 
x)^(1/2))^(1/2))/a^(5/2)+x*(a+c/x+b*(d/x)^(1/2))^(1/2)/a-3/2*b*d*(a+c/x+b* 
(d/x)^(1/2))^(1/2)/a^2/(d/x)^(1/2)
 
3.31.63.2 Mathematica [A] (verified)

Time = 0.79 (sec) , antiderivative size = 174, normalized size of antiderivative = 1.29 \[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=\frac {\sqrt {a} d \left (2 a-3 b \sqrt {\frac {d}{x}}\right ) \left (c+\left (a+b \sqrt {\frac {d}{x}}\right ) x\right )+\sqrt {d} \left (4 a c-3 b^2 d\right ) \sqrt {\frac {d \left (c+\left (a+b \sqrt {\frac {d}{x}}\right ) x\right )}{x}} \text {arctanh}\left (\frac {\sqrt {c} \sqrt {\frac {d}{x}}-\sqrt {\frac {d \left (c+a x+b \sqrt {\frac {d}{x}} x\right )}{x}}}{\sqrt {a} \sqrt {d}}\right )}{2 a^{5/2} d \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \]

input
Integrate[1/Sqrt[a + b*Sqrt[d/x] + c/x],x]
 
output
(Sqrt[a]*d*(2*a - 3*b*Sqrt[d/x])*(c + (a + b*Sqrt[d/x])*x) + Sqrt[d]*(4*a* 
c - 3*b^2*d)*Sqrt[(d*(c + (a + b*Sqrt[d/x])*x))/x]*ArcTanh[(Sqrt[c]*Sqrt[d 
/x] - Sqrt[(d*(c + a*x + b*Sqrt[d/x]*x))/x])/(Sqrt[a]*Sqrt[d])])/(2*a^(5/2 
)*d*Sqrt[a + b*Sqrt[d/x] + c/x])
 
3.31.63.3 Rubi [A] (warning: unable to verify)

Time = 0.35 (sec) , antiderivative size = 149, normalized size of antiderivative = 1.10, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.318, Rules used = {2065, 1693, 1167, 27, 1228, 1154, 219}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx\)

\(\Big \downarrow \) 2065

\(\displaystyle -d \int \frac {x^2}{d^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}}d\frac {d}{x}\)

\(\Big \downarrow \) 1693

\(\displaystyle -2 d \int \frac {x^3}{d^3 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}d\sqrt {\frac {d}{x}}\)

\(\Big \downarrow \) 1167

\(\displaystyle -2 d \left (-\frac {\int \frac {\left (2 \sqrt {\frac {d}{x}} c+3 b d\right ) x^2}{2 d^3 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}d\sqrt {\frac {d}{x}}}{2 a}-\frac {x^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{2 a d^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle -2 d \left (-\frac {\int \frac {\left (2 \sqrt {\frac {d}{x}} c+3 b d\right ) x^2}{d^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}d\sqrt {\frac {d}{x}}}{4 a d}-\frac {x^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{2 a d^2}\right )\)

\(\Big \downarrow \) 1228

\(\displaystyle -2 d \left (-\frac {\frac {\left (4 a c-3 b^2 d\right ) \int \frac {x}{d \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}d\sqrt {\frac {d}{x}}}{2 a}-\frac {3 b x \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{a}}{4 a d}-\frac {x^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{2 a d^2}\right )\)

\(\Big \downarrow \) 1154

\(\displaystyle -2 d \left (-\frac {-\frac {\left (4 a c-3 b^2 d\right ) \int \frac {1}{4 a-\frac {d^2}{x^2}}d\frac {2 a+b \sqrt {\frac {d}{x}}}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}}{a}-\frac {3 b x \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{a}}{4 a d}-\frac {x^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{2 a d^2}\right )\)

\(\Big \downarrow \) 219

\(\displaystyle -2 d \left (-\frac {-\frac {\left (4 a c-3 b^2 d\right ) \text {arctanh}\left (\frac {2 a+b \sqrt {\frac {d}{x}}}{2 \sqrt {a} \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}\right )}{2 a^{3/2}}-\frac {3 b x \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{a}}{4 a d}-\frac {x^2 \sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c d}{x^2}}}{2 a d^2}\right )\)

input
Int[1/Sqrt[a + b*Sqrt[d/x] + c/x],x]
 
output
-2*d*(-1/2*(Sqrt[a + b*Sqrt[d/x] + (c*d)/x^2]*x^2)/(a*d^2) - ((-3*b*Sqrt[a 
 + b*Sqrt[d/x] + (c*d)/x^2]*x)/a - ((4*a*c - 3*b^2*d)*ArcTanh[(2*a + b*Sqr 
t[d/x])/(2*Sqrt[a]*Sqrt[a + b*Sqrt[d/x] + (c*d)/x^2])])/(2*a^(3/2)))/(4*a* 
d))
 

3.31.63.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 1154
Int[1/(((d_.) + (e_.)*(x_))*Sqrt[(a_.) + (b_.)*(x_) + (c_.)*(x_)^2]), x_Sym 
bol] :> Simp[-2   Subst[Int[1/(4*c*d^2 - 4*b*d*e + 4*a*e^2 - x^2), x], x, ( 
2*a*e - b*d - (2*c*d - b*e)*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a, b, c 
, d, e}, x]
 

rule 1167
Int[((d_.) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_S 
ymbol] :> Simp[e*(d + e*x)^(m + 1)*((a + b*x + c*x^2)^(p + 1)/((m + 1)*(c*d 
^2 - b*d*e + a*e^2))), x] + Simp[1/((m + 1)*(c*d^2 - b*d*e + a*e^2))   Int[ 
(d + e*x)^(m + 1)*Simp[c*d*(m + 1) - b*e*(m + p + 2) - c*e*(m + 2*p + 3)*x, 
 x]*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, m, p}, x] && NeQ[m 
, -1] && ((LtQ[m, -1] && IntQuadraticQ[a, b, c, d, e, m, p, x]) || (SumSimp 
lerQ[m, 1] && IntegerQ[p]) || ILtQ[Simplify[m + 2*p + 3], 0])
 

rule 1228
Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c 
_.)*(x_)^2)^(p_.), x_Symbol] :> Simp[(-(e*f - d*g))*(d + e*x)^(m + 1)*((a + 
 b*x + c*x^2)^(p + 1)/(2*(p + 1)*(c*d^2 - b*d*e + a*e^2))), x] - Simp[(b*(e 
*f + d*g) - 2*(c*d*f + a*e*g))/(2*(c*d^2 - b*d*e + a*e^2))   Int[(d + e*x)^ 
(m + 1)*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x 
] && EqQ[Simplify[m + 2*p + 3], 0]
 

rule 1693
Int[(x_)^(m_.)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_.), x_Symbol 
] :> Simp[1/n   Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*x + c*x^2)^p, 
x], x, x^n], x] /; FreeQ[{a, b, c, m, n, p}, x] && EqQ[n2, 2*n] && IntegerQ 
[Simplify[(m + 1)/n]]
 

rule 2065
Int[((a_.) + (b_.)*((d_.)/(x_))^(n_) + (c_.)*(x_)^(n2_.))^(p_), x_Symbol] : 
> Simp[-d   Subst[Int[(a + b*x^n + (c/d^(2*n))*x^(2*n))^p/x^2, x], x, d/x], 
 x] /; FreeQ[{a, b, c, d, n, p}, x] && EqQ[n2, -2*n] && IntegerQ[2*n]
 
3.31.63.4 Maple [A] (verified)

Time = 0.28 (sec) , antiderivative size = 213, normalized size of antiderivative = 1.58

method result size
default \(\frac {\sqrt {\frac {b \sqrt {\frac {d}{x}}\, x +a x +c}{x}}\, \sqrt {x}\, \left (4 a^{\frac {5}{2}} \sqrt {b \sqrt {\frac {d}{x}}\, x +a x +c}\, \sqrt {x}-6 a^{\frac {3}{2}} \sqrt {b \sqrt {\frac {d}{x}}\, x +a x +c}\, \sqrt {\frac {d}{x}}\, \sqrt {x}\, b +3 \ln \left (\frac {\sqrt {\frac {d}{x}}\, \sqrt {x}\, b +2 \sqrt {b \sqrt {\frac {d}{x}}\, x +a x +c}\, \sqrt {a}+2 a \sqrt {x}}{2 \sqrt {a}}\right ) d a \,b^{2}-4 \ln \left (\frac {\sqrt {\frac {d}{x}}\, \sqrt {x}\, b +2 \sqrt {b \sqrt {\frac {d}{x}}\, x +a x +c}\, \sqrt {a}+2 a \sqrt {x}}{2 \sqrt {a}}\right ) a^{2} c \right )}{4 \sqrt {b \sqrt {\frac {d}{x}}\, x +a x +c}\, a^{\frac {7}{2}}}\) \(213\)

input
int(1/(a+c/x+b*(d/x)^(1/2))^(1/2),x,method=_RETURNVERBOSE)
 
output
1/4*((b*(d/x)^(1/2)*x+a*x+c)/x)^(1/2)*x^(1/2)*(4*a^(5/2)*(b*(d/x)^(1/2)*x+ 
a*x+c)^(1/2)*x^(1/2)-6*a^(3/2)*(b*(d/x)^(1/2)*x+a*x+c)^(1/2)*(d/x)^(1/2)*x 
^(1/2)*b+3*ln(1/2*((d/x)^(1/2)*x^(1/2)*b+2*(b*(d/x)^(1/2)*x+a*x+c)^(1/2)*a 
^(1/2)+2*a*x^(1/2))/a^(1/2))*d*a*b^2-4*ln(1/2*((d/x)^(1/2)*x^(1/2)*b+2*(b* 
(d/x)^(1/2)*x+a*x+c)^(1/2)*a^(1/2)+2*a*x^(1/2))/a^(1/2))*a^2*c)/(b*(d/x)^( 
1/2)*x+a*x+c)^(1/2)/a^(7/2)
 
3.31.63.5 Fricas [F(-1)]

Timed out. \[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=\text {Timed out} \]

input
integrate(1/(a+c/x+b*(d/x)^(1/2))^(1/2),x, algorithm="fricas")
 
output
Timed out
 
3.31.63.6 Sympy [F]

\[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=\int \frac {1}{\sqrt {a + b \sqrt {\frac {d}{x}} + \frac {c}{x}}}\, dx \]

input
integrate(1/(a+c/x+b*(d/x)**(1/2))**(1/2),x)
 
output
Integral(1/sqrt(a + b*sqrt(d/x) + c/x), x)
 
3.31.63.7 Maxima [F]

\[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=\int { \frac {1}{\sqrt {b \sqrt {\frac {d}{x}} + a + \frac {c}{x}}} \,d x } \]

input
integrate(1/(a+c/x+b*(d/x)^(1/2))^(1/2),x, algorithm="maxima")
 
output
integrate(1/sqrt(b*sqrt(d/x) + a + c/x), x)
 
3.31.63.8 Giac [A] (verification not implemented)

Time = 0.39 (sec) , antiderivative size = 218, normalized size of antiderivative = 1.61 \[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=-\frac {2 \, \sqrt {a d^{2} x + \sqrt {d x} b d^{2} + c d^{2}} {\left (\frac {3 \, b}{a^{2}} - \frac {2 \, \sqrt {d x}}{a d}\right )} + \frac {{\left (3 \, b^{2} d^{2} - 4 \, a c d\right )} \log \left ({\left | -b d^{2} - 2 \, \sqrt {a d} {\left (\sqrt {a d} \sqrt {d x} - \sqrt {a d^{2} x + \sqrt {d x} b d^{2} + c d^{2}}\right )} \right |}\right )}{\sqrt {a d} a^{2}} - \frac {3 \, b^{2} d^{2} \log \left ({\left | -b d^{2} + 2 \, \sqrt {c d^{2}} \sqrt {a d} \right |}\right ) - 4 \, a c d \log \left ({\left | -b d^{2} + 2 \, \sqrt {c d^{2}} \sqrt {a d} \right |}\right ) + 6 \, \sqrt {c d^{2}} \sqrt {a d} b}{\sqrt {a d} a^{2}}}{4 \, \sqrt {d} \mathrm {sgn}\left (x\right )} \]

input
integrate(1/(a+c/x+b*(d/x)^(1/2))^(1/2),x, algorithm="giac")
 
output
-1/4*(2*sqrt(a*d^2*x + sqrt(d*x)*b*d^2 + c*d^2)*(3*b/a^2 - 2*sqrt(d*x)/(a* 
d)) + (3*b^2*d^2 - 4*a*c*d)*log(abs(-b*d^2 - 2*sqrt(a*d)*(sqrt(a*d)*sqrt(d 
*x) - sqrt(a*d^2*x + sqrt(d*x)*b*d^2 + c*d^2))))/(sqrt(a*d)*a^2) - (3*b^2* 
d^2*log(abs(-b*d^2 + 2*sqrt(c*d^2)*sqrt(a*d))) - 4*a*c*d*log(abs(-b*d^2 + 
2*sqrt(c*d^2)*sqrt(a*d))) + 6*sqrt(c*d^2)*sqrt(a*d)*b)/(sqrt(a*d)*a^2))/(s 
qrt(d)*sgn(x))
 
3.31.63.9 Mupad [F(-1)]

Timed out. \[ \int \frac {1}{\sqrt {a+b \sqrt {\frac {d}{x}}+\frac {c}{x}}} \, dx=\int \frac {1}{\sqrt {a+\frac {c}{x}+b\,\sqrt {\frac {d}{x}}}} \,d x \]

input
int(1/(a + c/x + b*(d/x)^(1/2))^(1/2),x)
 
output
int(1/(a + c/x + b*(d/x)^(1/2))^(1/2), x)